Mon. Not. R. Astron. Soc. 000, [TH24] (2008) Printed 7 May 2008 (MN WI^ style file v2.2) 



Two-fluid models of superfluid neutron star cores 

N. Chamer 

Institut d'Astronomie et d'Astrophysique, Universite Libre de Bruxelles, CP226, Boulevard du Triomphe, B-1050 Brussels, Belgium 
7 May 2008 



ABSTRACT 

Both relativistic and non-relativistic two-fluid models of neutron star cores are con- 
structed, using the constrained variational formalism developed by Brandon Carter 
and co-workers. We consider a mixture of superfluid neutrons and superconducting 
protons at zero temperature, taking into account mutual entrainment effects. Lep- 
tons, which affect the interior composition of the neutron star and contribute to the 
pressure, are also included. We provide the analytic expression of the Lagrangian 
density of the system, the so-called master function, from which the dynamical equa- 
tions can be obtained. All the microscopic parameters of the models are calculated 
consistently using the non-relativistic nuclear energy density functional theory. For 
comparison, we have also considered relativistic mean field models. The correspon- 
dence between relativistic and non-relativistic hydrodynamical models is discussed in 
the framework of the recently developed 4D covariant formalism of Newtonian multi- 
fluid hydrodynamics. We have shown that entrainment effects can be interpreted in 
terms of dynamical effective masses that are larger in the relativistic case than in the 
Newtonian case. With the nuclear models considered in this work, we have found that 
the neutron relativistic effective mass is even greater than the bare neutron mass in 
the liquid core of neutron stars. 

Key words: stars: neutron - dense matter - hydrodynamics - relativity - equation 
of state 



INTRODUCTION 



The recent discovery of quasi periodic oscillations (QPOs) in the X-ray flux of giant flares from soft-gamma repeaters (SGR) 
may well be the flrst direct observational e vidence of neutron star oscillations. QPOs have been detected during the 2004 
December 27 giant flare from SGR 1806-20 Jlsrael et al.ll2005l:lstrohmaver fc Wattslbood : IWatts fc Strohmave jEooel 'l . during 
the 1998 August 27 giant fl are from SGR 1900-1-14 i Strohmaver fc Wattjboosi ) and during the 1979 March 5 event in SGR 



0526-66 jSarat et alj|l983l ). Those QPOs are usually interpreted as glo bal seismic vibrations trigg ered by magnetic crust 
quakes (for a recent review of these so-called magnetars, see for instance IWoods fc Thompsonll2006l and references therein). 
If this interpretation is conflrmed, the analysis of those Q POs can potentially reveal the interior composition of neutron 
stars, thus putting constraints on the theory of dense matter i Samuelsson fc Anderssonl 2007 '). Apart from QPOs in SGR, the 
rapid development of the g ravitational wave astronomy opens very exciting perspectives of directly observing neutron star 
oscillations in a near future ( Andersson fc KokkotaslboOSi ). In particular, accreting neutron stars in Low Mass X-Ray Binaries 
are expected to be detectabl e by the Advanced LIGO detectoi0 which is planned to be operational in a few years (see for 
instance the recent analysis of Watts et al.ll200^ . However the interpretation of these observations requires not only a detailed 
theoretical understanding of the dynamics of neutron stars, but also a consistent description of the different layers. 

A neutron star is mainly composed of three distinct regions: an outer cru st, an inner crust characterized by the presence of a 
neutron ocean and a liquid core which might be solid in the deepest regions ( Haensel. Potekhin fc Yakovlevll2006|'). Microscopi c 
calculations of dense nuclear matter suggest that the matter inside neutron stars is superfluid (IDean fc Hiorth-Jensenll2003l ). 
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This theoretical prediction is strongly supported by the observations of pulsar glitches i Bavm et al. 19691 : Anderson fc Itoh 
I975I ). Other in dications in favour of super fluidity, while less convincing, are provided by observations of neutron star thermal 



X-ray emission (jYakovlev fc Pethickll2004l ). One of the remarkable consequences of superfluidity is the possibility of having 
several dynamically distinct components. The electrically charged particles inside neutron stars are locked to gether by the 
interior magnetic field and co-rotate on very long time scales of the order of the age of the star I Easson 19791 ). The charged 
particles are rotating at the observed angular velocity of the star due to the coupling with the radiating magnetosphere 
and thus follow the long-term spinning-down of the star caused by the electromagnetic radiation. In contrast the neutrons 
being electrically uncharged and superfluid can rotate at a different rate. This naturally leads to considering the interior of a 
neutron star as a two-fluid mixture. As a result of the strong interactions between neutrons and protons, the two fluids are 
not completely independent but are coupled via mutual entrainment effects. These non-dissipative effects are known to affect 
significantly the frequencies of superfluid o scillation modes for which the superfluid neutrons and the charged particles are 
counter moving (lAndersson fc Comer 2001 ). T wo-fluid models of superflu i d neutron s tar cores including entrainment e ffects, 
have been proposed by Comer fc Jovnt (|2003l ). Carter. Chamel fc Haensel l|2005l , l2006l ) (see also lchamel fc Carterlliooj ) have 
shown how to describe in a unifled way, both the liquid core and the inner crust with in this two-fluid picture. The description 
of the outer crust requires a different treatment ( Carter. Chachoua fc Chamell 20 06|'). Eventually the differ ent layers have to 
be matched with appropriate boundary conditions at the interfaces ijAnderssonet al. 20021 : Lin et al. 2007 1. 

In this work, we have constructed relativistic hydrodynamical models of cold superfluid neutron star cores, calculating 
all the necessary micr oscopic coefficients with the same underlying microscopic model. The present work differs from that of 



Comer fc Jovnt (|2003h by improving the microphysics description of dense nuclear matter and by making the link between 



non-relativistic and relativistic models. In the first Section, we briefly review the convective variational formalism of multi- 
fiuid systems. This approach is employed to construct a two-fluid model of neutron star core in the Newtonian framework in 
Section [3l It is then shown in Section |4] how to generalize this model to relativistic fluids. In Section O we discuss the effects 
of entrainment in terms of dynamical effective masses. The consequences of superfluidity at the hydrodynamical scale are 
discussed in Section|6l Section[7]is devoted to the conditions of "chemical" equilibrium and to the composition of neutron star 
cores. In Section [S] the microscopic parameters of the two- fluid model are evaluated, using Skyrme effective nucleon-nucleon 
interactions. As an example, numerical results are shown in Section|9]for three particular Skyrme forces: the popular SLy4 force 
and the parameter sets LNS and NRAPR which were entirely constructed from realistic quantum many body calcu l ations . 
In Section [TOl results are compared to those obtained by the relativistic mean fleld theory applied by Comer fc Jovnt ( 20031 ) . 
We have cons tructed new relativisti c mean field models that yield a much better agreement with nuclear data than those 
considered by I Comer fc Jovnt (|2003f ). 



2 VARIATIONAL FORMALISM OF MULTI-FLUID HYDRODYNAMICS 



In the following, we will use Greek letters for the space-time indices . . . and Latin letters i,j, . . . for space indices. We 
introduce capital Latin letters x,y for distinguishing the various constituents and among them we will adopt the symbols i 
for leptons and q for nucleons (when several indices of the same species will be needed, we will add primes on the label as for 
instance q, q' , q" , etc.). We will apply the Einstein summation convention (i.e. repeated indices are summed) for space-time 
indices but not for the constituent labels. 

Let us co nsider a n arb itrary number of fluids that are interac ting with each other. We follow the variational forma lism 
developed by Carter ( 19891 ). which has been recently reviewed by Gourgoulhon ( 2006h and Andersson fc Comer (2007). In 
this approach, the basic fluid variables are the particle 4-currents of each fluid. The equation governing the dynamical 
evolution of each fluid is obtained from an action principle by considering variations of the fluid particle trajectories. Given 
a so-called master function A, which is the Lagrangian density of the system, and a set of 4-force densities acting on each 
fluid, the hydrodynamic equations take the very simple form 



where the vorticity 2-form cc/Ji is defined as the exterior derivative of the 4-momentum covector 

^ _ dA 

namely 



2Vr 



(1) 



(2) 



(3) 



It is understood that the partial derivative in Eq. ([JJ is taken with all other 4-currents being kept constant. For a strict 
application of the variational principle, the forces should separately vanish /5 — 0, which entails by contracting Eq. ([1} with 



the corresponding 4-current n^, that 



0. In this case, Eq. ([T| therefore reduce to 
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= 0. (4) 

Note that despite the vanishing of the individual forces , the fluids are not independent of each other in general and the set 
of Eqs Q are therefore not simply Euler equations. The couplings between the various fluids are hindered in the momenta 

The stress-energy tensor of the fluids can be expressed as 

X 

where ^ can be interpreted as a generalized pressure and is defined by 

X 

Let us emphasize that so far we have made no assumption with respect to the space-time geometry so that the above covariant 
expressions, based only on the exterior calculus, are valid both in (special and general) relativity and in the Newtonian limit. 
In the following sections we will show how to construct the Lagrangian density A in each case. 



3 NON-RELATIVISTIC TWO-FLUID MODELS OF NEUTRON STAR CORE 



We consider a uniform mixture with four constituents: neutrons, protons, electrons and possibly muons. Such a composition 
is expected to be found in the interior of low-mass neutron stars and in the outer core of massive neutron stars at densities 
above the crust-core transition density pcc ~ po/2 and below < 2 — 3po, where po — 2.8 x 10^'* g.cm~^ is the saturation density 
of infinite symmetric n uclear matter. Given the current unc ertainties on the composition of neutron star core (for a recent 
review see for instance Haensel. Potekhin fc Yakovlev 2006h . it is sometimes assumed for simplicity that this composition 
remains the same at higher densities. 

At densities below ~ 3 po, the nucleons are essentially non-relativistic. For instance, the sound velocity for the realistic 
model A18 + Sv + UIX* of Akmal et al. ( 19981 ). becomes comparable to the speed of light at densities around ~ 5po- For 
this model, such densities are reached in very massive neutron stars with a mass larger tha n 2M(7). Howev er the most 
precisely measured neutron s tar masses (in neutron star binaries) lie below 1.5Mq ILattimer fc Prakash 200?! '). Besides as 
shown bv iGlendenning (chapter 3, Section 4), the effects of General relativity are negligible at the microscopic scale. 



Moreover at the macrocoscopic scale, the fluid velocities are small compared to the speed of light. Indeed the velocity at the 
equator of the most rapidly spinning neutron stars is only about ~ 20% of the speed of light (taking 1 ms for the period and 10 
km for the radius). For the purpose of matching the microscopic nuclear model to the macroscopic hydrodynamical model, it 
will therefore be convenient to start with a local analysis considering non-relativistic fluids in the Newtonian framework. Such 
non-relativistic models can be also very us eful by themselves for studying qualitatively the dynamics of superfluid mixtures 
in neutron stars (jAndersson fc Comeill200lf ). In order to facilitate the correspondence be t ween r elativi stic and non-relativistic 
models, we will use the fully 4D covariant formalism developed by Carter fc Chamel 1 2004 . 2005al lb[). We will then show 
how to construct fully relativistic fluid models in Section |4] (as required for a General Relativistic description of the star). 
For simplicity we only consider the possible presence of the magnetic field by supposing that leptons and protons are co- 
moving, as discussed in Section [T] We therefore consider only two independent fluids: the neutron superfluid and the fluid of 
charged particles (protons, electrons and possibly muons). This two-fluid model includes the limit of non-super fluid neutron 



star cores since in this case, all the particles are essentially co-moving and can thus be treated as a single fluid (|Bavm et al 



19691) . Including the magnetic field is in principle straightforward. It has been recently shown that under some circumstances 
the magnetic field can even be va riationally taken in to account in the purely Newtonian context despite the non-Galilean 



invariance of Maxwell's equations ( Carter et al. 2006h . However taking into account the magnetic field as a dynamical field, 



implies a better understanding of the lepton dynamics as well as the proton superconductivity which is beyond the scope of 
the present model. 



Let us now briefiy review the 4-dimensional geometric structure of the Newtonian space-time (see ICarter fc Chamel 



20041 for a detailed discussion). Newtonian theory postulates the existence of a universal time t, leading to a foliation of the 
space-time into 3-dimensional hyp ersurf aces. Each of these spatial sections are fiat and are endowed with the 3-dimensional 
Euclidean metric, giving rise to the symmetric contravariant tensors rj'^'^ and rj^v by pushforward and by pull-back respectively. 
These tensors are not metric tensors since they are degenerate 



Tj^^U = , rypi.e" = 

where ti, — d^t and the "ether" fiow vector , normalized by the condition 



(7) 



(8) 
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characterizes a particular Aristotelian frame corresponding to the usual kind of 3+1 space-time decomposition. The particle 
4-currents introduced in Section [2] are given by = n^u^, where n-^ are the corresponding particle number densities, and 
the 4-velocities uj^ are defined by 

< = <+e^ t,< = 0. (9) 

with V-^ being the corresponding push forward of the usual 3-velocities in the given Aristotelian frame. It is easily seen from 
Eq. (IHJ that the 4-velocities are normalized as 

iM< = 1 • (10) 

We will neglect the small mass difference between neutrons and protons and write simply m for the nucleon mass, which 
we take equal to the atomic mass unit. For consistency we thus take the electron mass m'^ = 0. We also neglect the muon 
mass except at the microscopic scale for calculating the internal energy density. In the limit of small currents, the Lagrangian 
density A can be decomposed into a dynamical part Adyn, which depends on the particle currents n^, and a static part Ains 
which depends only on the particle densities given by 

= "-x^M • (11) 

The dynamical Lagrangian density (neglecting the contribution of the leptons) can be written as a quadratic form of the 
nucleon currents 

Adyn=iE^M.^'''«'- (12) 

9,9' 

The coefficients of the symmetric mobility matrix {q,q' = n,p for neutrons, protons respectively) are functions of the 
nucleon densities and will be evaluated in Section [H] Due to the Galilean invariance, the matrix elements are related to each 
other by 

Y^ng,K:'"'' =m. (13) 

9' 

Taking the partial derivative of the above equation with respect to the nucleon density leads to the following identity 

q',q" 

Since the left-hand side of Eq. p4|) has to be negative for any neutron and proton densities, we obtain the following inequalities 

driqi ' \ dUq J dUq dUq 

The non-diagonal coefficient IC"^ = IC" accounts for the non-dissipative entrainment effects and arises from the strong 
interactions between nucleons. If the nucleons could be treated as ideal Fermi gases, this coefficient would simply vanish 
K."" = 0. 

The static contribution Ains is related to the static internal (non gravitational) energy density Uh-is by 

Ains = —Uina — f/pot , (16) 

where p is the total baryon mass density 

p — mni, (17) 

and C/pot = pfj} is the gravitational potential energy density, (j) being the gravitational scalar potential. As discussed previously 
nucleons are not relativistic neither at the microscopic scale nor at the macroscopic scale. Nevertheless, in order to prepare 
the generalization to relativistic fluids in the next section, the nucleon rest mass energy density is included in t/ins. It can be 
shown that in the Newtonian framework this extra term has no effect o n the variational pr inciple (thereby on the dynamical 
equations of the fluids) provided the total mass current is conserved ( Carter et al. 20061 ). Let us point out that the total 
internal (non gravitational) energy density Uint is not simply given by f/ins, but contains an additional entrainment term 
defined by Ucnt = Udyn — C/i^in, where Udyn = Adyn and 



^X^*?^" r"'^'?""-'?'' (IS) 



is the total kinetic energy density. We thus have Uint ~ Uins + Ucnt- 

In summary, the non-relativistic Lagrangian density of the fluids is therefore given by 

A = Adyn + Ains , (19) 
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where Adyn and Aina are given by Eqs. ([12} and (|16|l respectively. The dynamical equations of the gravitational field can be 
obtained from the variational principle by simply adding the contribution 

Agrf = -^r;'"'(V^<^)(V.<^) (20) 

G being the gravitational constant. Considering variations of the total Lagrangian density Atot = A + Ag^f with respect to 
leads to the usual Poisson's equation 

r?''"V^V.(;!> = 47rGp. (21) 

Given the above Lagrangian density (|19|l . the momentum of each constituent is obtained from Eq. ([2}. The nucleon 
momentum is given by 



K = 5Z "•9' + *M ( ^ X] ^""^^ — n^,ng„ - - m<;/> j , (22) 

q' \ q',q" j 

where the chemical potential defined by 

Since the Lagrangian density depends only on the lepton densities = n'^tii {£ = e,fj. for electrons, muons respectively), the 
momenta of the leptons are time-like and are given by 

^1 = -i^Mf • (24) 
If the particles are all co- moving with the 4- velocity u'^, the nucleon 4- momenta take the familiar expression 

TT^ = mr]^„u'^ — t^ (^^"^^^ + ^9 + "^0^ (25) 

where v'^ — rj^^u'^u^ , using the identities (I13p and (|14p . 

The usual 3-momentum covector, denoted by 11-^, is defined by the Aristotelian spatial components of the 4- momentum 
covector n^, 

= n'i^l , (26) 

where rji^ is the space projection tensor defined by 

< ^Si;^ e^t. , (27) 

using the Kronecker unit tensor 5^. It follows immediately from Eq. (|8]) that eTI^ — 0. It is readily seen that the lepton 
3-momentum vanishes fl^^j — 0, while the nucleon 3-momentum is determined solely by the mobility matrix 

q' 

From the nucleon and lepton momenta Eqs. (|22p and (|24p respectively, we can obtain the generalized pressure 'i' of the 
fluids according to Eq. ((6)1. The kinetic part of the Lagrangian density Akin ~ f^kin does not contribute to the pressure 
which can thus be written as 

* = A.„.-^n,-^, (29) 

X ^ 

where Aint = A — Akin- Due to entrainment eff'ects, this internal Lagrangian density Aint of the fluids is not simply equal to 
the opposite of the total internal energy density (including gravitational contribution) Uint + Upot but is given by 

Ai„t = -Uint - Upot + 2f/cnt . (30) 

The gravitational potential energy density does not contribute to the pressure. As a result, the generalized pressure 'I' is the 
sum of a static term ^Pins given by 

*i--E"x^-f/in., (31) 



and a dynamical term given by 

dUon 



*e.t=-E<S^ + t^-- (32) 
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In the single fluid case, the static pressure 'Pins reduces to the ordinary pressure usuaUy denoted by P. Note that in muhi-fluid 
systems, if the static internal energy density is of the form 



(33) 



the static pressure can then be written as the sum of the partial pressures 



(34) 



with 



dm 



(35) 



In the general case however such a decomposition is not possible. The additional pressure term 'I'cnt vanishes whenever either 
the two fluids are co-moving or the fluids are non-interacting so that the non-diagonal coefficients of the mobility matrix 
vanish (no entrainment). 

Going back to the two-fluid model of neutron star cores, the static and entrainment contributions to the general pressure 
are given explicitly by 



Wins = n„ — h rip 



dn„ 



dn„ 



+ Ue 



and 



dn, 
2 drn 



n„ 



dn„ 



+ n: 



9 dml 



^ drip 



+ in„(m" — m) + ^np^m^ — m) 



(36) 



(37) 



respectively. Given the general pressure 'I' — ^'cnt -I- 'I'ins and the 4-momenta tt^ of the various constituents, the stress-energy 
tensor of the fluids can easily be obtained from Eq. ([5j . 



4 FROM NON-RELATIVISTIC TO RELATIVISTIC TWO-FLUID MODELS 

In the previous section, we have taken into account of the effects of the gravitational field on the fluids by including the term 
—p4> in the static internal Lagrangian density Ains. Alternatively as was first shown by Elie Cartan, the effects of gravitation 
can be taken into account in the structure of the Newtonian space-time itself, thus facilitating the comparison with General 
Relativity (see for instance chapter 12 of Misner et al. 1973h . This can be achieved by replacing the tensor r;^^ in Eq. (|12p 



for the Newton-Cartan space-time metric j^i, defined by 

■Ji^L, = Vi^f — 20t^t^ . (38) 

It is readily verifled, remembering Eq. (|13p , that this is indeed completely equivalent to adding the term —p4> to the Lagrangian 
density A. However using this approach, it becomes clear how to make the correspondence with General Relativity (thereby 
Special Relativity as well). For weak gravitational fields, the Riemannian metric g^i, of the relativistic space-time can be 
locally approximated by 

~ ri^„ — {(? -I- 2(j))t^t„ = 7,i„ — c^t^t„ . (39) 
This suggests to define the "dynamical" contribution to the relativistic Lagrangian density as 

Adyn = i ^ /C' {g^^nl^rig, + c^nqrig,) , (40) 

where rig — riqUq . The particle densities rix are now defined by 

= \J~at^i^ri^ri^/c, (41) 
with the 4-velocities normalized as 

5m.« = . (42) 

Note that Eq. 1)41^ with the normalisation (I42|l are consistent with Eqs. (|9} and Eq. (llOf) in the non-relativistic limit. Likewise 
with the definition (|40p . the expression (|12p is indeed recovered in the Newtonian limit. 

The relativistic expression of the corresponding "static" part is readily obtained by simply substituting the particle 
densities of the constituents in the Newtonian ether frame for the densities (|4ip in the rest frame of the corresponding 
particles. Using the identity (|13p and adding the internal energy density, the total relativistic Lagrangian density of the fluids 
can be expressed as 
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A = i g^^n^Tiq, + ^n\,mc - Uina , (43) 

1,1' 

where n^, — rin + Up is the baryon density. If the fluids are described in General Relativity, the Lagrangian (|43|l has to be 
complemented with the Einstein-Ifilbert contribution 

A.. = ^i?, (44) 

where R is the Ricci scalar associated with the metric g^^^. Variations of the action with respect to the metric (which involve 
not only variatio ns of the total L agrangian density but also variations of the space-time measure) leads to Einstein's equations 



(see for instance I Andersson fc Comer ,2007, ). 



The relativistic momenta of the nucleons can be expressed in a form similar to Eq. (|26|) as 

The non-diagonal components of the symmetric relativistic mobility matrix IC'' , are equal to those of the non-relativistic 
matrix IC'' , while the diagonal elements are given by 

j^qq ^ j^qq „ + J!l T^— — (cn„np + Qp^n^n") . (46) 

Uq C^Uq C^Uq OUq ^ ^ ' 

The relativistic momenta of the leptons take a very simple form 

= ^"^p 9t^^ ■ (47) 
If the constituents are all co-moving with the 4- velocity m^, the 4- momenta reduce to 

<' = ^"''ffA.- (48) 

With the momenta specified, we can obtain the generalized pressure ^ from Eq. (IS}. As for the non-relativistic case, ^ 
can be decomposed into an ordinary "static" part given by Eq. and an extra contribution ^tent due to entrainment which 
can be expressed as 

*ent = -E<^ + f^-. (49) 

where the entrainment energy density is now defined by 

C/cnt = ^E/C««'5M-n,^n;,. (50) 

9.9' 

Before concluding this section, let us remark that the relativistic Lagrangian density can be written in the very concise 
form 

A = Ao + Ai(x^-np), (51) 
where the coefficients Aq and Ai are given by 

Ao = -C/ins , Ai = -c'/C"^ , (52) 
and adopting the following notations 

X C ~ gfj,i/Tl^Tlp , (^3) 

n^c^ = -gp^^n^Un , (54) 
p^c = — g^^Tip . (55) 
Equat ion (|51|) is consistent with the expansion of the Lagrangian density in powers of (x^—np), suggested by Andersson. Comer fc Langloij 



( 2002ll . \t can be clearly seen on Eq. (|5ip . that in the absence of entrainment (i.e. Ai = 0) or in the case of co- moving fluids. 



the Lagrangian density reduces to the opposite of the internal energy density. 



5 DYNAMICAL EFFECTIVE MASSES 
5.1 Non-relativistic case 

If the nucleons were not interacting with each other, the mobility matrix introduced in Section [3] would be diagonal and we 
would simply have A^"" — m/rin and lO'^ = m/rip. Of course we know that nucleons are strongly interacting. This means 
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that the matrix IC'' does not have such a simple structure. It is convenient to define neutron and proton dynamical effective 
masses by 

m" = Tin IC"" , = Up IC^ (56) 

respectively. The deviations of ml from the bare baryon mass m therefore arise entirely from the nucleon-nucleon interactions. 
Let us point out that these effective masses depend on the nucleon densities and therefore vary with depth inside the neutron 
star. As a result of Eq. (|13[) . the non-diagonal coefficients of the mobility matrix can be expressed as 

n p 

jQ^p JQ^^ ^ ^* ^ ^* (5*^) 

Tip Tin 

It can be shown with stability arguments I Chamel fc Haensefcood ^ that the effective masses are bounded from below ml/m > 
Tiq/rii,, where Wb = n„ + rip or equivalently 

rih 

Since in neutron star core the effective masses are typically smaller than the bare nucleon mass (therefore IC"^ > 0), this 
inequality provides an upper bound for the largest possible strength of entrainment effects between the two fluids. The physical 
meaning of the dynamical effective masses defined by Eq (|56p . becomes clear when writing the expressions of the nucleon 
3-momentum covectors (1261) 



= m"vn,^ + {m — m")vp^ , (59) 

^1 = m^Vp ^ + {m — m^)v„ v . (60) 

where i/ = rj^pV^. This shows that the 3-momentum and the 3- velocity of a given nucleon species are not aligned whenever 
the dynamical effective masses differ from the bare nucleon mass, or equivalently whenever the non-diagonal coefficients of 
the mobility matrix do not vanish. 



5.2 Relativistic case 

By analogy with the definition (|56p . let us introduce relativistic nucleon dynamical effective masses by 

ihl = TlgiC"" , (61) 

where is the relativistic generalisation of the non-relativistic mobility matrix . Using Eq. (I46p together with l|56[l. we 
find 

= 1 2 ^1 2 ^ \^ '^"'^P + dpTrinrip ) . (62) 

m m mc mc' oriq 

with Hq the chemical potential defined by Eq. (|23p . It is easily checked that ml — > ml in the Newtonian limit. With these 
definitions, the neutron and proton 4-momenta can be explicitly written as 

tt'/j = [m"un + {Tn - m")up] g^^ , (63) 

= [m^Up + (m - m^)Un] g^i^ ■ (64) 

Note that the entrainment contributions involve the non-relativ istic effective masses (|56p . Eq. (|62|l is the generalization to 
interacting multi-fiuid systems of the effective mass introduced by Carter ( 19891 ) in the perfect fiuid case. Indeed in the absence 



of entrainment, K."^ = so that ml = m while the relativistic effective masses are given by 

^ = (65) 



This equation is identical to Eq.(1.66) in the lectures notes of Carter ( 19891 ). The physical origin of the difference between 



fhl and ml is that in relativity all forms of energy contribute to the mass. A remarkable consequence is that even massless 
particles can have a non-vanishing relativistic effective mass. This is for instance the case of leptons for which we have assumed 
m^ = 0. However Eq. (|47p show that the dynamical effective lepton mass is not zero but is given by 

m^, = — . (66) 

Note that even if leptons were not interacting (which we will actually suppose in Section |8] in order to evaluate the master 
function A), they would still have an non-zero effective mass due to the Pauli exclusion principle which prevents all the 
particles from occupying the lowest energy state with zero momentum (the chemical potential is then given by the Fermi 
energy of the lepton species I). 

For unbound nuclear systems like the liquid core of neutron stars, we have ^q > mc^. Besides if we assume that the 
strength of entrainment effects decreases with increasing density, i.e. 
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^<o, (67) 

OUq 

(this is actually the case for the models considered in this work, see Eq. (|109|l l. and using the Cauchy-Schwartz inequality 

- gti,^n,'inp < \/-g^,,nlin,'(^-gf,vn^n^ (68) 

it is easily shown that the relativistic effective masses are always larger than the non-relativistic ones. With the notations 
introduced at the end of Section[4l the relativistic dynamical effective masses can be expressed solely in terms of the parameters 
Ao and Ai of the relativistic Lagrangian density as 

= (n, - n,)^ - + l|^(np - x^) . (69) 

Note that the first order expansion of the relativistic Lagrangian density, Eq. (|51|l in powers of {x^ — np) leads to dynamical 
effective masses with a first order contribution proportional to (x^ — np). This velocity-dependent term vanishes in the 
Newtonian limit so that the non-relativistic effective masses are independent of the velocities to first order. The reason is that 
in relativity the particle number densities involve all the components of the corresponding 4-currents according to Eq. (|41[) 
while in the Newtonian case, the densities only depend on the time component of the 4-current through Eq. 



6 SUPERFLUIDITY 

In the previous sections, we have accounted for the superfluidity in neutron star core by assuming that two independent fluid 
motions could co-exist. However strictly speaking this assumption only requires perfect fluidity, i.e. the absence of viscosity 
and dissipative drag effects which damp the development of relative motions between the constituents. The distinguishing 
feature of a superfluid compared to a perfect fluid is the fact that it can be described by a macroscopic quantum wave function. 
This entails that in a superfluid the momentum circulation is quantized according to the Bohr-Sommerfeld quantization rule 

TT ^dx" = NTvh , (70) 

where h is the Dirac-Planck constant and is an integer, which simply follows from the requirement that the length of any 
closed path must be an integral multiple of the de Broglie wavelength of the condens ate formed of bound neutron pairs. This 
condition implies the existence of neutron quantized vortex lines in neutron stars ( Ginzburg fc Kirzhnits 19651 ). Assuming 



that the neutron vortices are arranged on a regular triangular array, the inter vortex spacing is given by 



where fin is the angular velocity of the neutron superfluid, which is approximately equal to the observed angular velocity 57 
of the star. 

In regions devoid of vortices, the neutron momentum circulation is equal to zero which implies that the neutron momentum 
can be written as 

= ^V,^" , (72) 

where the factor of 1/2 accounts for the fermionic nature of the neutrons and if" is the scalar phase of the condensate. 
Consequently, the corresponding vorticity 2-form locally vanishes 

zu^^ = . (73) 

However at length scales much larger than for which we are interested here, a fluid element is threaded by many vortex 
lines. Consequently the vorticity 2-forms do not have to vanish at this scale. Nevertheless, the superfluidity condition requires 
the existence of an average 4-velocity vector uJJ of the vortex lines such that the Lie derivative of the vorticity 2-form along 
u{5 vanishes 

u^£m;^,^0. (74) 
The above condition is satisfled if 

■"JSro^Ti. = . (75) 

Vortices are co-moving with the superfluid, unless forces act on them. Indeed it can be seen that Eq. (|75p with u!i = ul^ is 
consistent with Euler Eq. (|4]| obtained for the case /" — 0. It should be stressed that the above conditions (|73|l and (|75l) 
apply for either relati vistic or non-r elativistic superfluid. The presence of vortices can be explicitly included in the variational 
principle as shown bv lCartei] (|2000l ) and will not be further discussed here. 
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7 COMPOSITION OF NEUTRON STAR CORE 



The composition of the neutron star core is determined by the rates of transfusion processes which convert particles of different 
species into each other. While the baryon number is always conserved 



0, 



(76) 



where n'^ = n'^ ~\- Up, neutrons may be transformed into protons and vice et versa via electroweak processes. The fastest 



process is the direct Urea process 



n + ve, 



(77) 



where £ is electron or muon. When the beta equilibrium is reached, the two reactions occur at the same rate. In degenerate 
dense matter, this process is all owed for sufficiently large proton fractions owing to the requirement that both momentum 



and energy has to be conserved i Lattimer et al.|[l99lh . When these reactions are forbidden, the slower modified Urea process 
prevails 

n + N -* p+ +e + ue, 



(78) 



involving an additional spectator nucleon (neutron or proton). The relaxation time of these beta processes for npe matter, 

neglecting nucleon superfluidity, is approximately given by r'^' ~ 20T^'^ s an d r^*^-* ^ ' 

Urea processes respectively, where Tg is the temperature in units of 10^ K ( Yakovlev et al 
transformed into each other via the lepton modifled Urea processes 



months for the direct and modified 
Electrons and muons are 



20011 ) 



ll~ +X 



+ x- 



Ve + Vu 



(79) 



where X is either a nucleon or a lepton (the dire ct process is kin cmatically forbidden). The relaxation time associated with 
electromagnetic processes, of the order of 10~^^ s ( Easson &: Peth ick 1979), is much smaller than the characteristic time-scales 
of the neutron star phenomena considered here, so that the matter can be treated as electrically neutral. This condition reads 

Up ^ Ue + . (80) 

In the newly-born proto-neutron stars, the temperatures are of the order of ~ 10^^ K or higher so that the equilibrium 
is reached in a few microseconds for the modified Urea process or ten times less for the direct Urea. As the star cools down 
to temperatures ~ lO'' K after 10"^ — 10* years, the relaxation times rise dramatically to about 20 seconds for the direct 
Urea and several months for the modified Urea. Besides whe n the temperature falls b elow the critical threshold for the onset 
of superfluidity, the relaxation times increase exponentially ( Villain &: Haense] 2005 ). As a consequence, for the short time- 
scales ~ 1 — 100 milliseconds relevant for oscillations of mature neutron stars (like the recently observed QPOs in SGR), the 
composition of the star remains essentially frozen and the constituents can therefore be assumed to be separately conserved 



0, 



0, 



(81) 



which entails by Eqs. ([76} and (|80|) . that the other currents are also conserved V^nJ^ = and V^njJ — (remembering 
that the leptons are co- moving with the protons). Let us remark that the lepton number is not conserved unlike the baryon 
number, because the neutron star matter is transparent to neutrinos (except for the flrst few seconds after its birth into a hot 
proto-neutron star). 

The initial equilibrium composition of the neutron star core is obtained from the conditio n of electro neutralitv (|8pil and 
the conditions that the chemical affinities corresponding to the above processes sh ould vanish ( Carter fc Chamel 2005b|). The 
chemical affinity of a given reaction S is defined by IjCarter fc Chamelll2005bl ) 



(82) 



wh ere and £^ a r e the relevant particle creation numbers and the energies per particle, respectively. As pointed out 
bv lCarter fc Chamell (|2005bh . the problem arises of determining the reference frame with respect to which the energies £^ 
have to be measured when some of the constituents (here neutrons and charged particles) are moving with different velocities. 
Since the relative velocity between the two fluids is expected to be small compared to the fluid velocities, we assume for 
simplicity in this section that the particles are all co-moving with 4- velocity u'^. It is then natural to deflne the energy per 
particle as f ^ = — u'^tt-^,. In the Newtonian ease, using Eq. (|25p we thus have 



(83) 



where the chemical potential of a particle species x is deflned by Eq. (|23|) . Since the mass is conserved in any chemical reaction 
H involving non-relattvistic particles, the corresponding chemical affinity (|82p reduces to 



(84) 
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In the relativistic case, the energy per particle obtained from Eq. (148 |l is given by 

£^ = Mx - (85) 



so that Eq. (|84p is valid for both relativistic and non-relativistic particles. 

We assume that neutrinos have escaped from the star so that we set f^*^ = 0. Both the direct and modified Urea processes, 
respectively (|77|) and (|78|) . have the same affinity given by 



_^Urca^^„_^p_^e^ (86) 

The electron-muon transfusion reactions (|79|l are characterised by the affinity 

A'"'=£^-£''. (87) 

The composition of the core at a given baryon density nb, can then be determined by solving the equations ^^'^'^^ = = A'^'^ 
under the constraint (|80|) . Using (|84|) . this leads to 

= Me + , (88) 

Me = Mm • (89) 

If the matter is in equilibrium, the static pressure depends only on the baryon density rib and can be written in the 
concise form (valid in both the relativistic case and the Newtonian limit) 

*ins = ribMn - f/ins (90) 

where fin is the neutron chemical potential evaluated at the equilibnum neutron and proton densities, associated with the 

baryon density nb. Let us emphasize that Eq. (|90p is only valid for neutron star matter in equilibnum. In the general case, 
the static pressure is given by Eq. (|31|l . 



8 EVALUATION OF THE MICROSCOPIC PARAMETERS 

The internal static energy density can be decomposed into several contributions 

f/ins {n„, rip. Tie, 71^} = C/n Tip} + (Jcoui {rip, He, u^} + [/l {n^} + (7l {n^} , (91) 

where C/n is the nucleon part, C/coui is the Coulomb part and U-l the lepton (kinetic) part. The Coulomb energy arises from 
lepton-lepton, lepton-proton and proton-proton interactions. It is of purely quantum origin since the classical contribution 
coming from Poisson's equation vanishes as a result of electro neutrality. We neglect the lepton-lepton interactions (but not 
the lepton-proton interactions) and we approximate the proton Coulomb energy by the Hartree-Fock exchange energy of 
non-relativistic point-like charged particles 

f/coui {np} = -|e^ . (92) 

Unlike nucleons, leptons are relativistic at the microscopic scale (note however that this does not imply that their collective 
motion is relativistic at the macroscopic scale of the fluid description) . Their kinetic energy is thus given by that of an ideal 
relativistic Fermi gas {£ — e,jj,) 

he 



C/l {ne} 



xt(2xj + l)y^2;| + 1 - In ^a::^ + ^ x'j + 1^ 



(93) 



where \i — fi/mic is the Compton wave length and the dimensionless parameter xi is defined in terms of the Fermi wave 
number 

fcpf = (37r^nf)^/^ (94) 
by 

xi = Af fcpf . (95) 
Since the electron mass is set to zero, the electron energy density is obtained by taking the limit Xe —* +oo of Eq. (|93p yielding 



^^K}^^, (96) 
where ppe = ftfcpe is the electron Fermi momentum. 
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The strong interactions among nucleons are d escribed by an effective Hamiltonian with a two-body force of the Skyrme 
type (|Bender et alJbooi Istone fc ReinhardlboOTi ) 



v{ri,r2} = to{l + xoPa)S{r} + —{1+xiP^ 

is, 



k'^S{r} + 5{r}k^ 



+ -^{l+x-sP^)S{r}rn,{Rr + iWo{ai+a2) ■ k' 



t2{l+X2Pa)k' ■ 5{r}k 
X 5{r}k 



(97) 
(98) 



where r = ri — 1*2, R = {ri +r2)/2, cti and 0^2 are Pauli spin matrices, k — — i(Vi — V2)/2 is the relative wave vector, fc' is the 
complex conjugate of k acting on the left, and Pa = (I+cti •CT2)/2 is the spin-exchange operator. The first term represents the 
attractive part of the nucleon-nucleon interaction. The next two momentum dependent terms are associated with the finite 
range of the interaction. The density dependent term proportional to ts corresponds to the strongly repulsive short range part 
of the interaction and simulates the effects of three body interactions (Vautherin & Brink 1972). The l ast term which leads 



to spi n-orbit coupling in finite nuclei does not contribute in uniform matter. In principle, as shown bv lNeeele fc Vautherin 



I 1972h . ;his effective interaction can be derived from the "bare" nucleon-nucleon interaction by expanding the nucleon density 
matrix in relative and centre of mass coordinates, r, R respectively. In practice however, the parameters are usually determined 
by fitting experimental data and/or results of microscopic many-body calculations in infinite uniform nuclear matter using 
the bare nucleon-nucleon interactions. Such kind of zero range effective forces are valid whenever the inter particle spacing 
is much larger than the range of the nuclear interactions. This condition is satisfied at densities below ~ 3po- Nevertheless 
since these effective forces are usually constrained to reproduce the high density equation of state of nuclear matter, it is not 
completely unreasonable to apply them at densities above ~ 3po- Finite range effects of the nucleon-nucleon interaction as 
well as relativistic corr ections are somehow t aken into account phenomenologically by the fitting procedure. For instance, the 
parametrizations SLy (IChabanat et al.lll997l ) have been specifically constructed for neutron star studies by fitting a "realistic" 
equation of state of neutron m atter up to very high densities ~ lOpo- Besides, it is worth mentioning that soon after such 
effective forces were introduced. ICameronI (1953) applied them to neutron stars and showed that the maximum mass ~ 2Mq 
is compatible with the scenario of neutron star formation from supernova explosions. The main limitation of these effective 
forces is th at they describe only nucl eonic degrees of freedom. At densities above ~ 3po, other particles like hyperons are likely 
to appear jHaensel. Potekhin fc Yako vlcv 200^). Nevertheless, let us remark that for the most precisely measured neutron 
star masses in binary ra dio pulsars, their cen tral densities lie below ^ 3 — 4po depending on the equation of state (see for 
instance chapter 6 from Haensel et al. ( 20061 ) and in particular their figure 6.3). Our main motivation for using such kind of 
effective Hamiltonian is the perspective of a unified treatment of the interior of neutron stars, including not only the liquid 
core but also the solid crust whose microscopic description starting from the bare nucleon-nucleon interaction is not feasible. 

Assuming that the mat ter in neutron star co res is not polarized for the densities of interest as suggested by many-body 
calculations (see for instance Bombaci et ahlbood ) , the nucleon energy densi ty associated with the force (|97p can be calculated 
using the method outlined in the classic paper of I Vautherin fc Brink! 1 19721 ) and is given by an expression of the form 



Un {n„,np} = nt,mc? + - — Tb -f -Bin^ + B2{n\ + nl) + B^n^r^ + BA{nn,T, 
2m 

where 



+ UpTp) + Bsn^^" + Ben^ (n: 



2 I 2n 



|(37^^)^/^4/^ 



(99) 



(100) 



are respectively the neutron and proton kinetic energy densities (in units of Ir^ /2m), and — t„ + Tp. The B-coefficients are 
related to the parameters of the force by the following expressions 

to /-^ _|_ ^0' 



Bi = 



R - *° 
B2--J 



B4 = - 



i2 1 + -f 



tl { Xl + - 



t2 



1 

X2 + - 



Bs = 



^fl + ^i 
12 V 2 



Bk = 2:3 H — 

12 V 2 



(101) 
(102) 
(103) 
(104) 
(105) 
(106) 



Analytic expressions of the effective masses m" and for the force (|97p have been recently obtained by Chamel fc Haense] 

1 20061 ). ' ntroducing the parameter /Ss = 2mBz/ti? , the coefficients of the mobility matrix, given by Eqs. (|56[l and (|57[). can be 
expressed as 
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^„„^ml+^ (107) 
n„ 1+ PaUh 

Tip 1 + /JsWb 

J^np ^ j^pn ^ jn—^ . (109) 

Note that as a consequence of the isospin symmetry of the nucleon-nucleon interactions, we have /C""'{n„, n^} = IQ"^ {np,nn} 
and K,"^ (therefore Ai) does not depend on the matter composition but only on the total baryon density rib = n„ + rip. 
This means that entrainment effects are not affected by the various chemical reactions that may occur inside the core, as 
discussed in Section [T] In the high density limit pzUn 3> 1 and /Sarip ^ 1, all the elements of the mobility matrix become 
equal m/nh- As a consequence, the non-relativistic effective masses tend to ml/m —> Uq/nh. This asymptotic limit 

which corresponds to the strongest entrainment effects (see the discussion of Section ^ is never reached in neutron star core 
for the nucleon-nucleon interactions considered in this work, since m/(5z is typically of the order of ~ lOpo (see Table|3}. 



9 CHOICE OF EFFECTIVE MICROSCOPIC HAMILTONIAN 



We have selected effective forces according to the following criteria. First of all, the chosen forces have to yield reasonable 
values of the "semi-empirical" saturation properties of infinite uniform symmetric nuclear matter, namely the equilibrium or 
saturation density no (or the mass density po = nam), the binding energy per nucleon 

[/N{no/2,no/2} 2 



the symmetry energy coefficient 



_ 1 9 

" 2 dP 



with / = (ri„ 



9n, 



rib , 

I — Q^n^—riQ 

rip) /rib, and the incompressibility modulus 
a' /(7N{rrb/2,nb/2} 



rib 



(110) 



(111) 



(112) 



"b="0 



Global fits to essentially all the a vailable experimenta l nuclear mass data yield uq ~ 0.16 fm a„ ~ — 16 MeV, ~ 28 — 35 



MeV and ~ 220 - 240 MeV (|Lunnev et al.lbooal '). Due to the strong interactions, the mass of the individual nucleons in 
nuclear matter is different from the bare mass and can be written as 
m 



= (1 + /)- 



m 

mt 



= (1-/)- 



(113) 



in which ml and m* are the so-called isoscalar and isovector effective masses respectively (see for instance iFarine et al.ll200ll ). 
The isovector effective mass is a crucial microscopic input since it controls directly the strength of entrainment effects in 
neutron-proton mixtures. Indeed the parameter /Ja which determines the mobility matrix, Eqs. (|107p , (|108p and (|109p . is given 

by 



m 

m* 



1. 



(114) 



In principle, this isovector effective mass can be determined from measurements of the giant isovector electric dipole resonance 
in finite nuclei (consisting of relative motions between neutrons and protons). Nevert heless estimates are model dependent 
providing values mj/m ~ 0.7—1 at saturation density (see in particular the discussion of iLunnev et al ]|2003l in Sect. III-B-5-e). 
Microscopic many-body calculations in infinite uniform nuclear mat ter starting from the bare nucleon-nucleon interaction lead 
to an isovector effective mass around m*/m ~ 0.7 (see for instance Zuo et al. 20061 ). Besides we consider only those effective 
forces that have been constrained to fit the uniform infinite neutron matter equation of state. Otherwise these effective forces 
could not be reliably extrapolated to the neutron rich matter inside neutron star core. 

T he main deficiencies o f effective forces is the existence of insta bilities that are not found by microscopic calcula- 
tions ( Margueron et al. 2002 : Agrawal et al. 2004; Lesinski et al. 20061 ) . Especially many Skyrme forces predict a spurious 
ferromagnetic transition in neutron matter above some critical densities. We thus require that no such instabilities occur in 
the density range of interest p < 3po by imposing that the dimensionless Lan dau parameter, usually noted Go, be greater 
than —1 in neutron matter (following the analysis of Margueron et al. 20021 ). It turns out that this criterion is very re- 



strictive. Several forces that reproduce reasonably well both the saturation properties of symmetr ic nuclear matter and the 
neutron matter equation of state do not pass this test. For instance, the parametrization RATP ( Ravet et al. 19821 ). which 
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Table 1. Parameters of the chosen Skyrme forces. The units of energy and length are MeV and fm respectively. 





SLy4 


LNS 


NRAPR 


to 


-2488.91 


-2484.97 


-2719.7 


tl 


486.82 


266.735 


417.64 




-546.39 


-337.135 


-66.687 


ts 


13777.0 


14588.2 


15042.0 


a 


1/6 


1/6 


0.14416 


XQ 


0.834 


0.06277 


0.16154 


Xl 


-0.344 


0.65845 


-0.0047986 


X2 


-1 


-0.95382 


0.027170 


X3 


1.354 


-0.03413 


0.13611 



Table 2. B-cocfhcients of the chosen Skyrme forces. The units of energy and length are MeV and fm respectively. 





Bi 


B2 


B3 


B4 


B5 


Be 


SLy4 


-1763.39 


1660.1 


32.473 


49.3128 


1925.34 


-2128.55 


LNS 


-1281.48 


699.233 


44.5497 


-39.0001 


1194.94 


-566.35 


NRAPR 


-1469.69 


899.595 


85.0067 


-55.9836 


1338.81 


-797.364 



was the first attempt to construct an effective force for astrophysical applications, predicts that neutron matter becomes 
spin polarized slightly above saturation density ~ 0.175 fm"'' (the density for the onset of instability is obtained by solving 
Go = —1). Like wise the forces SkM and Skyrme 1', which have been applied to study dense matter in neutron stars and 
supernova cores ( Bonche fc Vautherin 19821 : Lattimer et al. 19851 : Lassaut et al. 1987 : Lorenz et al. 19931 ). yield a ferromag- 
netic transition density in neutron matter ~ 0.2 1 2 fm~^ and ~ 0.256 fm~^ respectively. We have found that only the forces 
of the Saclay-Lyon group (jChabanat et al.l Il997l . Il998bl lal) and the recent parametrization LNS (|Cao et al.l 120061 ) satisfy all 
the above conditions. They predict a ferromagnetic instability in neutron matter like the other forces, but at significantly 
higher densities ~ 3 — 4po which we do not consider in this work. The force LNS seems the most appropriate to describe 
neutron star core since it was constructed so as to reproduce recent results of microscopic diagrammatic calculations (based 
on Brueckner theory) of infinite uniform nuclear matter with two- as well as three-body forces. In particular, this effective 
force not only fits well the energy per nucleon in symmetric and asymmetric nuclear matter, but fits also the nucleon effective 
masses for different asymmetries and different densities which directly determine the entrainment coefficients as previously 
discussed. Nevertheless, the SLy forces, which were constrained to reproduce some properties of finite nuclei (apart from the 
other constraints that we imposed) would be preferable if not only the liquid core but also the crust layers would have to be 
described with the same underlying microscopic H amiltonian. Besides the equ ation of state of neutron star matter with the 
force SLy4 has been tabul ated and widely ap plied (jHaensel fc Potekhirj|2004l ). For comparison, we hav e also considered th e 
parametrization NRAPR I Steiner et al. 2003) since it was adjusted on the realistic equation of state of Akmal et al. ( 19981 ). 
Nevertheless this force leads to a ferromagnetic instability at rather low density p < 2po- 

The parameters of the forces and the associated B-coefficients introduced in Section [S] are given in Tables [1] and [2] 
respectively. The nuclear matter properties predicted by these forces are summarized in Table |3l Figure [T] shows the binding 
energy per particle in uniform infinite neutron matter defined by E/A = Uj^{nn,0}/n„ — m(? . Let us stress that the LNS 
force was fitted to the latest results of many body calculations with two- and three-body forces, while the forces of the Lyon 
group were adjusted to reproduce an older neutron matter equation of state based on variational methods. 

Figures|3l|4]and[5]show the equilibrium composition of cold neutron star matter, composed of neutrons, protons, electrons 
and muons, obtained by solving Eqs. (|80|) . (|88|l and (I89|l . The figures show the electron, muon and proton fractions, defined 
respectively by rie/rib, n^/n^, and Up/n^,, as a function of the mass-energy density p = Uins/(? , which is approximately given 
by p ~ nbm for rib < 3no. All forces predict the appearance of muons at rib — 0.12 fm^'^ or p ~ 2 x 10^* g.cm~^). The 
forces LNS and NRAPR yield similar composition. They both predict a slightly larger (resp. smaller) proton fraction than 
the force SLy4 above (resp. below) po- This can be understood by remarking that Eq. (|88|l can be approximately written 
as jMuther et al.lll987l ) 



\l/3 



45{nb}(l-2xp 



(115) 



where the symmetry energy 5 {rib} defined by 



5{nb} = 



j7N{nb,0} i7N{nb/2,nb/2} 



> 0, 



(116) 
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Figure 1. Binding energy per particle of uniform infinite neutron matter for the SLy4, LNS and NRAPR efi^ective forces. 



70| 




[fm ] 



Figure 2. Symmetry energy cS{nb} (in MeV) for the SLy4, LNS and NRAPR eff'ective forces as a function of the baryon density 
= n„ + rip. 



represents the cost in (nuclear) energy per particle to replace protons by neutrons in symmetric nuclear matter. From Eq. (|115p . 
we have 

As can be seen on Figure (2] the forces LNS and NRAPR yield a larger (resp. smaller) symmetry energy than the force SLy4 
above (resp. below) the saturation density no (note that the symmetry energy coefflcient Os ~ S{no}). 

Figure [S] shows the static pressure (|90|l of npe/i matter in equilibrium as a function of the mass-energy density p = i7ins/c^- 
In figures [T] [8] and [9l we compare the effective masses defined by Eq. ([56} for the two Skyrme forces. In both cases the neutron 
effective mass is close to the bare nucleon mass while the proton effective mass is significantly reduced. This is consistent with 
the inequality (|58p which implies that 



!^-^>7, (118) 
m m 

with / = {un — np)/nh. The relativistic effective masses defined by Eq. (jGlf) are shown on Figures [10] Illl and ll2l For simplicity 
we have considered that neutrons and protons are co-moving so that the effective masses are given by 
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Figure 3. Equilibrium fractions n~g^/n\^ of protons (p), electrons (e) and muons (^) inside neutron star core predicted by the SLy4 
effective force. 




Figure 4. Equilibrium fractions n^/nb of protons (p), electrons (e) and muons (/i) inside neutron star core predicted by the LNS 
effective force. 




Figure 5. Equilibrium fractions Uy^jn^^ of protons (p), electrons (e) and muons (/i) inside neutron star core predicted by the NRAPR 
effective force. 
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Table 3. Properties of infinite uniform symmetric nuclear matter for the chosen Skyrme forces, no (fm~^) is the nuclear saturation 
density, (MeV) is the binding energy per nucleon of infinite symmetric nuclear matter, as (MeV) the symmetry energy coeflicicnt, 
Kcx> (MeV) the compression modulus, m*/m, and m*/m are the isoscalar and isovector efi'ective masses respectively, Uf (fm~^) is the 
density at which a ferromagnetic instability occurs in neutron matter. Realistic values o f these parameters a re no ~ 0.16 fm-3, ay ~ -16 
MeV, as ~ 28 - 35 MeV, Koo ^ 220 - 240 MeV, m*/m ~ 0.6 - 0.9, m*/m ~ 0.7 - 1 l|Lunnev et al.ll2003h . 





no 


flu 


as 




m*slm, 


m*/m 


rif 


SLy4 
LNS 
NRAPR 


0.160 
0.175 
0.1606 


-15.97 
-15.32 
-15.86 


32.0 
33.4 
32.79 


229.9 
210.9 
225.7 


0.696 
0.827 
0.695 


0.801 
0.728 
0.605 


0.59 
0.62 
0.28 
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Figure 6. Static pressure of neutron star matter in equilibrium for the Skyrme effective forces and for the BL2 relativistic mean field 
model. 



These relativistic effective masses are significantly different compared to the non-relativistic ones. In particular both nuclear 
forces predict that tlie neutron relativistic effective mass is larger tfian tfie bare mass. Moreover, from Eqs. pi8|l and 
we have 



m 



mi 
m 



m 



— + > /- 
m mc' 



Me 



(120) 



so that the splitting of tlie relativistic effective masses is larger than that of the non-relativistic ones since /i^ > 0. 
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Figure 7. Neutron and proton effective masses, respectively mj/m and mj/m, defined by Eq. 15611 . in neutron star matter in equilibrium 
for the SLy4 effective force. 
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Figure 8. Neutron and proton effective masses, respectively m" /m and m^/m,, defined by Eq. I)56|l. in neutron star matter in equilibrium 
for the LNS effective force. 




Figure 9. Neutron and proton effective masses, respectively mj/m and m^/m, defined by Eq. I I56I I. in neutron star matter in equilibrium 
for the NRAPR effective force. 
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Figure 10. Relativistic neutron and proton effective masses, respectively mj/m and mj/m, defined by Eq. II61I I. in neutron star matter 
in equilibrium for the SLy4 effective force. Neutrons and protons are co-moving. 
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Figure 11. Rolativistic neutron and proton effective masses, respectively m,"/m, and m^/m, defined by Eq. in neutron star matter 

in equilibrium for the LNS effective force. Neutrons and protons are co-moving. 




Figure 12. Relativistic neutron and proton effective masses, respectively m"/m and mj/m, defined by Eq. 1611 1. in neutron star matter 
in equilibrium for the NRAPR effective force. Neutrons and protons are co-moving. 



10 COMPARISON WITH RELATIVISTIC MEAN FIELD MODELS 

A few years ago, IComer fc Jmr^ l|2003l ) developed relativistic two-fluid models of s uperfluid neutron st ar cores. They have 
determined the master function A using the eflFective relativistic mean field theory l|Glendenning||2000l ). In the model they 
considered, the interactions between nucleons arise from the exchange of two massive mesons: the scalar meson a with mass 
nicr and the vector meson tj^ with mass ttIcj. The former accounts for the long range attractive part of the nucleon-nucleon 
interaction while the latter gives rise to the short range repulsive part. In the nuclear field theory, particles are described by 
a microscopic Lagrangian density £ (not to be confused with the macroscopic Lagrangian density A of the fluids) given by 
(using units c = h = 1) 



C = V[7^i-Df,. - mD]i> + ^[(9''cr)(9^cr) - m^cr^] - ^uj^^lo^v + ^ml.Lo^uj^ 



(121) 



with the nucleon field tp and the antisymmetric tensor LUf_ii, = 9^0;^ — duUi^j, (7*^ denote the Dirac matrices and ip = ?/'^7'') 
The nucleon-meson couplings are introduced in the gauge covariant derivative 



Du 



da 



(122) 



and in the Dirac eflective nucleon mass 



(123) 
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Figure 13. Equilibrium fractions Uy^jni, of protons (p), electrons (e) and muons (^) inside neutron star core predicted by tlie BL2 
relativistic mean field model. 



where and are dimensionless coupling constants. The field equations, which actually only depend on the quantities 
Cct = Sa/n-a and = gt/rn^ in uniform infinite matter, are solved in the Hartree approximation (exchange terms are 
neglected) ignoring the contributions of antiparticles (the so-called no Dirac sea approximation). 

As for the non-relativistic energy density functional theory d iscussed in Section [9l the free parameters of the model 
have to be determin ed by fitting to some nuclear matter properties. IComer fc Jovnd (2003) adopted two parameter sets from 
Glendenninj (|200oh (chapter 4, table 4.4). However these parameters were not obtained for the Lagrangian density given by 



Eq. H12ip but for a more elaborate class of models which includes scalar self-interactions as well as the vector-isovector rho 
meson. As a result, b y dropping these extra terms in the Lagrangian density C, the resulting a — uj models considered by 
Comer fc Joto3 (|2003l ) predict unphysical nuclear matter properties, as can be seen in Table[5](see the discussion in Section[9]). 



Besides these models predict that neutron matter is bound as shown in Figure 1141 unlike quantum many body calculations 
using realistic nucleon-nucleon interactions. Note also that these models predict that neutrons and protons have the same 
(Dirac) effective mass mp for any nuclear asymmetry I = (n„ — np)/ni,, in contradiction to microscopic calculations (see for 
instance van Dalen fc Fu"ch3l2007j^ especially their fi gure 4). As such, these models are therefore unsuitable for applications to 
neutron stars as pointed out by lGlendennin3 (|2000h . In order to reproduce the properties of finite nuclei and infinite nuclear 
matter with the same level of accuracy as the non-relativistic effective forces discussed in Section [9l other mesons mus t be 
included. Besides non-linear self-meson interactions must be introduced (see for instance chapter 4 of iGlendenninjbood ). 

For the present time, we will restrict the discussion of the entrainment effects to the a — lu model. The models considered 
bv lComer fc Jovnd (2003) can be significantly improved (keeping in mind the inherent limitations of such models) by simply 
refitting the parameters. With only two free parameters Co- and Ct^ , only two of the symmetric infinite nuclear matter properties 
listed in Table|3]can be fitted exactly. We have constructed three new parameter sets BL1-BL3 by fixing the saturation density 
to no = 0.16 fm~'^ and (i) the binding energy per nucleon to — —16 MeV for BLl, (ii) the Dirac effective mass to m_D = 0.7m 
for BL2, (iii) the symmetry energy to — 28 MeV for BL3 (the fitting procedure did not converge for = 30 MeV). The 
parameters of these new models are given in Tabled As can be seen in comparing Table [S] and |31 the overall agreement with 
empirical nuclear data is still very poor refiecting the lack of flexibility of these a — uj models. The most important constraint 
for application to neutron stars is to reproduce at least the eq uation of state of ne utron matter. We have thus constructed the 
parameter set BL4 by fitting the realistic equation of state of Akmal et al. ( 19981 ). The result of the fit is shown in Figure [T4l 
as well as the predictions of the other mean field models. This should be compared with the results of non-relativistic effective 
forces in Figure[T] Note that the parameter sets GLI, GLII and BL3 predict incorrectly the existence of bound neutron matter. 
From these four models, it seems that the best compromise is achieved for the parameter set BL2, yielding reasonable values 
of the saturation density, compression modulus, Dirac efi'ective mass (which is an important quantity for entrainment effects 
as discussed in Section [9]) together with a fairly good fit of the neutron matter equation of state. Note however that this 
model is still very crude compared to the SLy4 or LNS Skyrme forces presented in Section [Q] In particular, the values of the 
symmetry energy as and the binding energy per nucle on ay, which are two b asic nuclear matter properties, are unrealistic. 

We have applied the general expressions derived by Comer fc Jovnt ( 20031 ) within the a — ui mean field models, to evaluate 
the entrainment parameters and to compare the results with those obtain ed using the non- r elativ istic effective energy density 
functional theory. We have greatly improved the models considered by Comer fc Jovnt ( 2003h (i) by refitting the meson 
coupling constants leading to a better agreement with nuclear data as discussed previously, and (ii) by including muons which 
affect the composition of neutron star core and contribute to the pressure. Leptons are treated as ideal relativistic Fermi gases 
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Figure 14. Energy per particl e of uniform infinite neutron matter for relativis tic a — uj mean field models. The cu rve labelled APR is the 
"realistic" equation of state of lAkmal et al. I lll998l '). using the ana lytical fit oflHeiselberg fc Hiorth-JensenI 1I2OOOI) . The rest mass energy 
has been subtracted out. GLI and GLII are the models used by I Comer fc JovntI | |2003D . BL1-BL4 are the parameter set constructed in 
this work. 




Figure 15. Relativistic neutron and proton effective masses, respectively mj/m and mj/m, defined by Eq. II6II I. in neutron star matter 
in equilibrium for the BL2 (solid lines) and GLI (dotted lines) relativistic mean field models. Neutrons and protons are co-moving. 



as discussed in Section |9] The parameters Ao and Ai in the expansion of the master function A, introduced in Section |4l are 
given by 

Ao = Ajo - C/l {rz^} , Ai = -^|o (124) 

where .4|o and A|o are given by Eqs.(63) and (A8) respectively in the paper of lComer fc JovntI (|2003h . Note that we have added 
the muon contribution in the Lagrangian density. Using the nucleon chemical potentials given by Eqs.(A9) and (AlO) of that 
paper, together with the lepton chemical potentials defined by Eq. (|23|) . we have determined the equilibrium composition of the 
neutron star core assuming co-moving particles as discussed in Section[7] Results are shown in Figure[T3]for the parameter set 
BL2. The proton fraction is very small at low densities unlike that predicted by non-relativistic effective forces. As discussed 
in Section [21 this can be understood from the very small (incorrect) value of the symmetry energy at saturation density 
(see Table [5]). As can be seen in Figure (6] the equation of state is however similar to that obtained for the non-relativistic 
effective nucleon-nucleon interactions. The reason is that matter in neutron star core is almost pure neutron matter and all 
models SLy4, LNS and BL2 reproduce reasonably well the neutron matter equation of state (see Figure HI and I14p . 
The relativistic effective masses introduced in Section U] are given by 

'rn" — nnB\o, fh^ — npC\o (125) 

where B\o and Cjo are given by Eq. (64) and (65) respectively in the paper of IComer fc JovntI l|2003l ). As shown in Figure [151 
the neutron effective mass predicted by the model GL2 is larger than the bare nucleon mass while the proton effective mass is 
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Table 4. Pa rameters of the rr — uj rel ativistic mean fi e ld mo dels discussed in this work. GLI and GLII are the models proposed by 
lGlendenniii3 (2000) and employed by IComer fc Jovnti (|20oi) after dropping scalar self-interactions and the rho meson. The models 
BL1-BL4 are new parameter sets introduced in this work (the coupling constants are given in fm^). Note that for the baryon mass we 
have taken the average of the neutron and proton masses as Comer & Joynt (2003). 







ci 


GLI 


12.684 


7.148 


GLII 


8.403 


4.233 


BLl 


14.6063 


11.0544 


BL2 


9.3353 


7.2624 


BL3 


23.2707 


14.4061 


BL4 


7.79346 


6.06748 



Table 5. Properties of infinite uniform symmetric nuclear matter for a — uj relativistic mean field models. The quantities shown are 
the sa me as those introdu ced in Table [3] Realistic values are no ~ 0.16 fm-3, a„ ~ -16 MeV, 28 - 35 MeV, ~ 220 - 240 

MeV llLunnev et al.ll2003l) . Note that recent relativistic many body calculations by Ivan Dalen fc Fuchsl 0003) predict that the Dirac 
effective mass at saturation is m^j ~ 0.7m. GLI and GLII are the models emploved by iGomer'fc'^ovnq (12003 ) while models BL1-BL4 
are new parameter sets introduced in this work. 





no 


a„ 


as 




mo/m 


GLI 


0.28 


-66.43 


34.69 


2117.5 


0.385 


GLII 


0.41 


-66.69 


41.29 


2252.0 


0.407 


BLl 


0.16 


-16 


20.09 


674.0 


0.543 


BL2 


0.16 


-1.70 


16.29 


249.0 


0.700 


BL3 


0.16 


-72.24 


28 


2217.4 


0.338 


BL4 


0.14 


1.7358 


13.88 


130.5 


0.774 



smaller, as obtained for non-relativistic models. For comparison we have also plotted the eflFective masses obtained with the 
model GLI considered by Comer fc Jovnt ( 20031 ) . The neutron effective mass obtained with this parameter set is decreased at 
low densities compared to the ordinary mass in contradiction to previous results. This is a consequence of the fact that the 
model GLI predicts (incorrectly) that neutron matter is bound, as shown in Figure [141 so that fin < mc^. 



11 CONCLUSION 



The recent detection of QPOs in SGR (most likely associated with seismic vibrations triggered by magnetic crust quakes) 
and future observations with gravitational wave detectors, offer new possibilities to probe the interior of neutron stars and to 
test the theories of dense matter. Nevertheless the reliable identification of the various oscillation modes calls for a consistent 
theoretical description of the star. As a first step towards this goal, we have constructed fully self-consistent relativistic two-fluid 
models of neutron star cores, composed of superfluid neutrons and a conglomerate of protons, electrons and possibly muons. 
The mutual entrainment effects between the two fluids, resulting from the strong nucleon-nucleon interactions, are properly 
taken into account. We have determined the expression of the Lagrangian density in the variational framework developed 
by Brandon Carter and co-workers. We have also shown how to make t he correspondence with non-relativistic models by 
applying the 4D covariant formulation of Newtonian hydrodynamics of ICarter fc Chamell (1200J) . We have determined all 
the coefficients of these models consistently with the same microscopic approach. In the perspective of describing not only 
the liquid core of neutron stars but also the crust layers, we have employed the nuclear energy density functional theory 
which has been already successfully applied to study both isolated nuclei and infinite nuclear matter. As an example, we 
have calculated the composition, the equation of state and the entrainment matrix of neutron star core for three different 
nuclear mo dels: the popular SLy4 mo del for which the equation of state of both the crust and the liquid core has been already 
tabulated (jHaensel fc Potekhinll2004l ) and the more recent LNS and NRAPR models which have been entirely constructed 
from recent realistic many bod y calculations. For com parison, we have also considered relativistic a — uj mean field models 
that have been first applied bv lComer fc Jovnt We have improved their models by refitting the parameters in order 

to obtain a better agreement with nuclear data, but still we could not reach the same level of accuracy as the non-relativistic 
models mentioned above. This would require the introduction of additional meson fields, especially the rho meson. Besides 
self-meson couplings should be taken into account. Numerical calculations with both effective energy density functionals and 
relativistic mean field models have shown that the dynamical effective nucleon masses arising from entrainment effects are 
smaller than the ordinary mass in the Newtonian case. Relativistic effects increase effective masses, since all forms of internal 
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energy contribute to the mass. A rather unexpected consequence which has not been usually discussed in the literature, is 
that relativistic effective masses can be even larger than the bare mass in the liquid core of neutron stars. 
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